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This paper describes a new method for measuring the duration of mode- 
locked picosecond laser pulses. It is similar to the Two Photon Fluorescence 
and Second Harmonic Generation methods in that it measures not the pulse 
duration directly but rather an autocorrelation function of the signal. It 
has the advantage in that it can be used for very low-power repetitive light 
signals. 

I. INTRODUCTION 

Two methods for determining the duration of picosecond pulses 
have been reported. The earlier method, which utilizes second har- 
monic generation (SHG) in nonlinear crystals, was reported independ- 
ently and almost simultaneously by M. Maier, W. Kaiser and J. A. 
Giordmaine; 1 J. A. Armstrong; 2 and H. P. Weber. 3 The second method, 
which utilizes the two-photon fluorescence (TPF) in certain dyes, was 
first reported by J. A. Giordmaine 4 and his co-workers. Measurements 
of pulse durations based on these methods have been reported by num- 
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erous authors. Extensive references to these papers are included in a 
review article by A. J. DeMaria, W. H. Glenn, Jr., M. J. Brienza and 
M. E. Mack. n 

Neither of these methods measures pulse duration directly. Both, 
in fact, measure the autocorrelation function G{t) of the intensity 
I(t), namely 



G(t) = f /(/)/(/ + r) dt. 

J — CO 



A careful measurement of the contrast ratio between G (t) at the peak 
of the pulse and G{r) between pulses is necessary in order to deter- 
mine whether or not the measured values represent the duration of a 
pure AM pulse. (This was first pointed out for TPF by H. P. Weber 
and independently by J. R. Klauder, M. A. Duguay, J. A. Giordmaine 
and S. L. Shapiro.) 7 If the contrast ratio is not the correct value for 
pure AM pulses, little concerning the signal waveform can be reliably 
inferred from a knowledge of G (r) . 

The beat frequency detection (BFD) method proposed in this pa- 
per measures a different quantity, namely 

H k (r) = m(r) + \F k (-r) 
where 



l'\(r) = I \ f*" &*(t)&(t + r) exp (-i2x| kj 



dt 



and 8(0 is the electric field of the optical signal, 1/T is the pulse repeti- 
tion rate and k/T is the frequency of the harmonic chosen for analysis. 
It too suffers from the fact that if the correct "contrast ratio"* is not 
observed, the result is ambiguous. 

If the proper "contrast ratio" is observed, then 8(0 is real (i.e., there 
is no angle modulation of the signal) and in cases of practical interest 
the pulse duration is short compared with the period T. Thus, for 
small k, kt/T <K 1 over the region where 8(0 contributes significantly 
to the integral and 



I'\(-r) = F t (r) = \r f" 8*(/)8(/+r) 

1 J-r/2 



dt 



Thus H k (T) is essentially the autocorrelation function of 8(0 when 
there is no angle modulation. It is therefore as good a measure of pulse 
duration as is the autocorrelation function of I(t). 

* This term will be defined in Section III. 
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Under certain conditions there is reason to believe that the pulses are 
"chirped." That is, they have the form 

8(0 = 8„(0 exp \j(u t + at 2 )} 

where S o (0 is a real, slowly varying envelope function. The measured 
quantity H k (r) has an interesting and distinctive behavior in this case. 
This is considered in Section III. It should be possible in some cases to 
determine by the BFD method how much, if any, "chirp" is present. 
Neither TPF nor SHG is capable of providing this information. 

II. DESCRIPTION AND ANALYSIS OF THE BEAT-FREQUENCY DETECTION 
METHOD OF PULSE-DURATION MEASUREMENT 

The output S(«) of a laser is a periodic function which can be written 
in the form 

8(0 = £ A n cos {(no> + u )t + <t>»} (1) 

n~— oo 

where 2tt/co = T is the period, co is the center frequency of the spectrum 
and A n is the amplitude of the nth line in the spectrum. If the phases <*>„ 
are related by 

4> n = n3> ( 2 ) 

where <S> is any constant, the signal 8(0 is a sequence of pulses with no 
angle modulation and the pulse envelope has the minimum possible 
duration for the set \A n }. Such a signal is called "mode locked". If 
equation (2) is not satisfied, angle modulation (and a wider pulse) 
occurs. In general, if 4> n is a set of random numbers a "noiselike" (but 
periodic) signal is obtained. The pulse duration of the signal 8(0 is 
therefore dependent on \<f>„\. 

Consider the circuit shown in Fig. 1. The laser output 8(0 is divided 
into equal parts by the 3 dB hybrid mirror (beam splitter) and the two 
parts are recombined in the same mirror after one path has suffered a 
delay t relative to the other. Beam transforming lenses may be required 
in order to assure alignment of the phase fronts of the two interfering 
beams. Mirror M 2 is mounted on a track so that the delay in one arm of 
the device can be varied over a range corresponding to one period of the 
laser signal. Mirror M\ is mounted on a piezoelectric acoustic modulator 
so that its position can be varied periodically over a distance of a few 
microns for reasons which will be explained below. The output of the 
hybrid is therefore 
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Fig. 1 — Diagram of the BFD apparatus. 
(8(/)+ S(/- r» 

00 

= £ A n [<xm {(nu + o )t + <l> n \ + cos {(nu + (o )(t — t) + *,}]. (3) 

n-— oo 

If this signal is now applied to a product demodulator (photodetector), 
one obtains as output: 

V(t) = [8(0 + 80 - r)] 2 



cos — + cos <\[n + - ]w + o) 



}■] 



(4) 



where 



Bj" = A n+K A n , 0™ = n+ « - <f> n 
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and terms in 2o> £ have been ignored as has a time shift of amount t/2. 
If this signal is now passed through a bandpass filter which is cen- 
tered at fcw and is sufficiently narrow that only the « = k term in the 
sum is passed, one obtains 



V k (t) = Z £»*' cos {hot + 8™ 



COS -y + COS Sin + 7JKJ t co rr 



(5) 



It should be noted that since d' n 0) = for all n, V (t) is independent 
of {<(>„} and independent of t. 

V (t) = E Bi m [l + cos {(ruo + a» )r}]. (6) 

One can readily understand this result by realizing that the V h {t) 
term in general results from the beats between all possible pairs of 
lines separated by an amount ha. In particular V {t) results from 
"beats" between each spectral line and itself. Since no beats between 
different lines contribute to V (t), no relative phases between lines 
influence the result. 
Returning to the general case, square law detection of V k (t) gives 

U k (r) = 2 | V k (r) | 2 

= Z Z B ( n k) Bl k) cos (9™ - O 

n m 

• |cos 2 —■ + h cos (n - m)coT + h cos { (n + m + /c)co + 2co } r 
+ cos -jp cos j(n + g jw + toorT 

+ cos M m + -Jco + o) >t • (') 

We observe that the entire w dependence is contained in the last three 
terms. If one of the mirrors in the interferometer is scanned through 
an optical wavelength, the contribution of these terms is cancelled out 
(since over an optical wavelength their average value is zero). This 
is accomplished by mounting one of the mirrors on a piezoelectric 
transducer and vibrating it. The resultant signal (after passing through 
a suitable low-pass filter) is given by 

ffiM = E T,B ( n k) Bi k) cos K fc > - O 



• cos 2 -y + \ cos (n — m)wT • 



(8) 



6 



THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1971 



We pause to evaluate equation (8) for two extreme cases. Consider 
the mode-locked case A, = n* which implies d'J c) - 6™ = 0. This gives 



U k (T) moda = cos" 

lucked 



k 



■cot 



l£#.]' + * I T, B n exp Qnur) 



(9) 



In this equation and in the remainder of this paper the superscript (k) 
on B n and B n is to be understood. The mode-locked case represents one 
extreme in the sense that it provides the shortest possible pulse consistent 
with a particular set of amplitudes \A n ). The "random case", i.e., 
<}> n a random variable uniformly distributed on [ — x, t], represents the 
opposite extreme — the expectation value of the signal power is constant 
in time. The expectation value of U k (r) for this case can be obtained by 
observing that the expectation value of cos {6„ — 8 m \ when the !<£„} 
are uniformly distributed on [0, 2t) is given by 

(cos (6„ - O) = 1 ^ & = (for all n, m), 
= 1 if n — m (for all k) , 
= otherwise. 



This gives: 



<ff.M) = 



' cos2 (¥) + i]? 5 "' 



k ?± 0. 



(10) 



Equations (9) and (10) are illustrated in Fig. 2. 

We now return to a consideration of the general case. Equation (8) 
can be written in the form 



U k (r) 



where 



cos" 



h 



2 ®n 

2 C3„ exp (inur) 



+ J 



CO 

^2 (P>„ exp (—jutcr) 



(11) 



Similarly set 



Then 



(B„ = R n exp (jd n ). 
a n = A n exp (j<j} n ). 



(B„ = a„ +t a: 
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(^B n ) : 




Fig. 2— Expected experimental results for the four cases discussed in the text. 

Consider the expression : 



X) (B, exp (j2r=i 



T 



where 



Define 



Then 



= I E a„ T *a*exp(j2^ 

I i r T/2 i t \ 

= \±J E*(-l)E(x - t) exp [fak J,) 



E(t) = 2 «n exp (j=^n 



dl 



F k (x) = I ~ j E*(t)E(t + .r) exp \-j%rk ^ 

f/,(r) = cos 2 (^)f,(0) + |F fc (r) + \F k {-r). 



(12) 
(13) 
(14) 
(15) 
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From equation (15) we can draw an interesting comparison be- 
tween the BFD method and the TPF and SHG methods. It is clear 
from equation (14) that F k (0) is the power spectrum of the intensity 
and therefore the autocorrelation function of the intensity is simply 



G(r) = J^FMrnpy-fak* 

Therefore the quantities measured by all three methods depend only 
on the quantities F k {r). The results of the BFD method depend on 
the values for a single k and all t; the results for TPF and SHG on 
the values for t = and all k. From equations (12) and (14), we have 



(16) 



X) (B„ exp (;'27r^n 



FM = 
while from equation (1) 

Kt) = 

Thus, /Mr) and I(t) are formally equivalent with the transformations 

B„ «-+ A n , 



2 «« exp [j2m^ 



Comparing equation (13) with (1) shows that 

&(t) = Re \E(t) exp (ju t)}. 

That is, E(t) is the (complex) signal from the laser with the carrier 
frequency omitted. 

For the mode-locked case, the pulses are very short compared with 
the repetition rate l/T. Thus for at least the first few k = 1, 2, • • • we 
can write (2*kt)/T « 1 for those values of t for which E(t) is significantly 
different from zero. This gives (since from equation (16) F(t) is even 
when the B n are all real) : 

ft(-r) - F„(r) = |; f n E*(t)E(t + r) dt 

I I J-T/2 

Thus the first term in equation (15) is just a raised cosine with period 
T/k, but the sum of the second and third terms is the square of the 
autocorrelation function of E(t). This sum provides a good measure 
of the pulse duration when the signal is mode locked. 

Since this technique, like TPF and SHG, possesses some ambiguity 
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if the signal is not mode locked, it is important to be able to distin- 
guish mode-locked behavior from non-mode-locked behavior. In order 
to do this, consider the quantity 

oo 2 



O k (0) = I 

3 

2 



n- — oo 

CO 



mode- locked case, 



O oo 

= - E Bl random-phase case. 

^ 71--00 



k*0 



For k = 0, n = and <B„ = B n is real for any set {<£„}. Thus 

3 



t/o(0) = 



ZB n 



for all cases. By considering the k = (dc) term the quantity 

CO 

n-— oo 

can be measured experimentally. Then if 

tf*(0) ^ tf (0) k = 1, 2, 3 • • • 
the signal is not mode locked. 

III. SPECIAL CASES 

It is instructive to consider how much U k (Q) differs from U o (0) for 
various cases. To this end we consider three special situations, namely, 
a pure FM signal, a model of partial mode locking, and a "chirped" 
pulse. Consider first the sinusoidally modulated pure FM signal: 

£(/) = cos {u t + <f>s'mcot} 

00 

= 2 J n (4>) cos {(wo + no))t}. 

n-— oo 

This is just a special case of equation (1) with 

A n = J n (0), 4> n = allw. 
The resulting expression for U k (j) is therefore 

U k {r) = cos 2 kj f [ n E JMJ n+ S)\ 

1 I °° 
+ o Z) J n (<l>)Jn +k (<l>) exp (jmor) 
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This can be reduced, by means of a well-known Bessel function 
identity to 



U k (r) = |/l(20 sin y)- 



Note that U (0) = \ whereas U k (0) = 0, k = 1, 2, 3, • • • . Also, the 
largest value that U k (r) k = 1, 2, • • • can assume is more than 9 dB 
below £?o(0) for the FM signal. 

The simplest model for partial mode locking is to let the phase #„ of 
the nth spectral line of the laser be uniformly distributed on the interval 
[—a, a] = [— eir, en]. Both a Monte Carlo calculation and an analytical 
calculation were made. The analytical calculation is described below. 
The results of both calculations are presented in Fig. 3. 

Consider equation (8). In order to determine the expectation value 
(ensemble average) (U k (r)) we must compute the values of (cos (0„, — 6„)). 
Three cases arise (for k 5* 0) : 

Case I: n — m 

(cos (0„ - O) = 1. 

Case II: n = m ±/c 

(cos (6 n - O) = (cos (<f> n+2k - 2</> n+k + 4>„)> m a. 

Case III: All subscripts distinct 

<cos (0 n - d m )) - (cos (0 n+t - n + <f> m+k - <f> m )) m b. 

Then 

U k (r) = b E Z £ B £ m [cos 2 (^f) + h cos (n - m)a>r] 

(*r)+fl ?** 

+ 2a £ fU?„ +t [cos 2 (^) + | cos (tor)] • 
This can be rewritten in the form 

U k (r) = ft COS 2 (^f ) 

+ G 2 + MS- 8 " cosn^r) 2 + (E J?„ sin fkor) 2 ] 



+ 



nf'm 



COS 
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Fig. 3(a)— (U(t)) and (Pit)) for ideal mode locking, (b)— (U(t)), (P(t)) and 
three representative curves from the Monte Carlo calculation of U(r) and P(t) 
for e = 5 (c) — (C'(t)>, (Pit)) and three representative curves from the Monte 
Carlo calculation of U(r) and Pit) for <= = 0.7. (d)— (U(.t)). (P(t)) and three 
representative curves from the Monte Carlo calculation of U(t) and P(t) for 



= 1. 
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where 

G. - KZB n ) 2 + (1 - 6) IB 1 . + 4(a - 6) £ BA +4 , 

n n n 

G* - 1(1 - b) Z# - (a - 6) ZB n B n + k . 

We now turn to the problem of evaluating a and b. It can be shown* 
that 



<cos (x)) = n 



sua €,• 



where a; = JZj.i t { , 

ti is uniformly distributed on [— e, , e,], 
{ ti } is a set of statistically independent random variables. 

From this, one immediately obtains : 

/sin a\ 2 /sin 2a\ /sin aY , v 

a = [—) \-*r) = {—) «»<«>■ 

6 = (smMX 



In particular 
(O k (0)) = G i + G 2 + $b\ZB n \ 2 

- 16 I Z A I 1 + 1(1 - b) ZBl + 3(a - 6) £BA + * . (17) 
Let 

R k = {vmm \T,B n \ 2 ) = ^^ 

(l-b)ZB 2 n Y,B n B n + k 

= 6+ (xbs +2(a ~ 6) (z^r- (18) 

n n 

Note that for a = (mode-locked case) R k = 1 while for a = it 
(random-phase case) 



(ITbJ 



^* — / V"* n \2 = Pk 



This result, even though it is not new, is derived in the Appendix. 
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R k can be thought of as the "contrast ratio" in analogy with the TPF 
case. Consideration of the above results shows that R k tt 1 for a < 0.5 
rad. If we assume p k = 0.1 (a reasonable value for about 10 spectral 
lines contributing to the signal) and neglect the third term in equation 
(18) [Note that this term is always negative— therefore its omission 
makes the result somewhat pessimistic] we obtain the following results 



a (rad) 


R k 





l 


0.5 


0.86 


1 


0.55 


1.5 


0.28 


2 


0.14 


2.5 


0.10 



Thus the "degradation" in "contrast ratio" becomes quite pronounced 
by the time a exceeds one radian. In this model the expectation value 
of the actual pulse power is given by 



^■1-7 



E Al + 



sin a 



^ A n exp (jnut) 



Thus as a departs from zero the expectation value of the power separates 
into two terms, one a constant term proportional to £ Al , the other a 
replica of the mode locked pulse diminished in amplitude by the factor 
sin 2 a/a. At a = 1 radian (where R k is beginning to fall off as sin 4 a/a 4 ) 
the pulse term in <P(0) is down only 1.5 dB and falling half as rapidly. 
The expectation values of Ui(t) and P(t) are plotted in Fig. 3 for 
the gaussian case 



[-©"] 



A. = exp |_-^ 

where / is the pulse repetition frequency and / is related to the band- 
width of the pulse. For this figure a value of 0.1 was used for the 
quantity ///<>. The results of a Monte Carlo type calculation are also 
illustrated in this figure. 

Finally we consider the case of "chirped" gaussian pulses, i.e., a 
repetitive train of pulses of the form 



EU) = exp 



i-m 



exp (jfif) . 
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This can be represented by the series 

E{t) = £ A n exp [j(no>t + yn*)] 

where A n = exp[— %(ria>/£l) 2 ]. The following result is derived for 
arbitrary A„, however. Here <f>„ = yn 2 and 9„ = y(n 2 + 2nk + k 2 — n 2 ) 
= y k(2n + k). 
Consider 



5„ exp 



2irtn \ 



W 






2 4 n 4„ +t exp [yyfc(2n + k)] exp (j 



2*j»\ 



7 1 / 



2 4„A n+ * exp ^ -y^j 

where rr = [^ + (T/Trjfcy]. 

Comparing this with equation (12), we see that this "chirping" results 
in a decrease in the amplitude of the raised cosine term and a splitting 
in the t term in equation (15) . In fact 

0M m w (ffyjg) + iFk (m + r ) + ^( r _ tm 

where in this result F k {x) is evaluated using the unchirped pulse since 
the phase term was handled explicitly. 

This result is readily understood when one recalls that ^(t) is 
formally equivalent to I(t) with this substitution 

A„ <-> B n , 

<f> n *-* B n . 
When A n is gaussian, B n is also, but when <f> n is quadratic in n, 8 n is 
linear. A phase shift which is linear in frequency corresponds to a 
shift of the origin in the "time" domain. This result is illustrated 
in Fig. 4. 

The curve for y = 0.05 corresponds to the case in which the pulse 
duration is ten times the reciprocal bandwidth of the signal. This is 
comparable to the observed values of pulse duration and bandwidth 
for Nd-Glass lasers. Thus this technique might provide a useful tool 
for investigating the amount of chirp on pulses from these lasers. 

It is of interest to consider the second-order statistics of U(t) in order 
to estimate the amount of fluctuation to be expected in the random 
non-mode-locked case. In order to derive the second-order statistics, 
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Fig. 4 — U(t) for a gaussian pulse with various amounts of chirp. 

it is necessary to compute [£/*(t)] 2 and to observe that 

(cos {(0 n - 6 m ) + (0„. - d m .)}) = 1 n = m, n' = m', 

= 1 n = m' , m = n' , 
= otherwise, 
(cos {($„ - O - (0»> - m -)}> = 1 n = m, n' = m' , 

= 1 n = n' , m = m', 
= otherwise. 
After some tedious manipulation one obtains the variance 



ES [cos 2 (*f ) + 



\ cos (n — m)cor p^5l . 



This result is plotted in Fig. 5 for the signal illustrated in Fig. 3a. 

It is also instructive to consider a somewhat different experiment. 
If we repeat the above experiment for the case where either the phase 
fronts of the signals do not coincide or the beams do not overlap at 
the photodetector we obtain, instead of equation (4) , 

V(t) = | 8(0 | 2 + | S(* - r) | 2 

= Z Z B n cos [hot + d n \ cos^p 
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Fig. 5 — Fractional rms deviation for the two statistical cases discussed in the text. 



and 

V k (t) = £ B n cos {hot + 9 n } cos ^ 

n ** 

which upon square law detection becomes 

U k (r) = 2 | y< T) | 2 - cos 2 yZI B n B m cos (0 n - m ). 

& n m 

This gives the following result. 

Uniformly distributed random phase: 



(U k (T)) = cos 5 



(if)?* 



4 ( ha', 



- oosM^JlE' E'BX). 



Mode-locked case <£„ = n<l») 

P,(r)-eo9 , (^)(EBJ'. 



This case is identical in form with the Hanbury Brown and Twiss 8 
experiment. 

The results for the mode-locked case and the random-phase case are 
summarized in Table I and illustrated in Fig. 2. 
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IV. EXPERIMENTAL VERIFICATION OF THE METHOD 

A scaled experiment was performed using a He-Ne laser with a c/(2t) 
frequency of 100 MHz. The pulse duration was measured by the BFD 
technique and compared with the actual pulse shape as observed by 
means of a high-speed photodiode and a sampling oscilloscope. Figure 6 
shows the results of the experiment and Fig. 7 shows the oscilloscope 
display of the pulse. The pulse duration of 0.7 ns as determined from 
the BFD experiment is in excellent agreement with the value measured 
with the oscilloscope. 

V. CONCLUSIONS 

The experiment proposed here should be capable of distinguishing 
between mode-locked and non-mode-locked behavior of the laser. If 
mode-locked behavior is observed, a measure of the pulse width is 
available from the shape of the dc voltage U(t) versus t curve. 

The experiment is relatively simple to instrument and should be 
readily applicable to any CW laser whose pulse repetition rate is 
sufficiently low to allow construction of narrowband "IF" circuits at 
that frequency. 



100 




-60 



-20 O 20 

ZS£ IN CENTIMETERS 



Fig. 6 — Experimental results for a mode-locked He-Ne laser. 
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Fig. 7 — Mode-locked pulse from He-Ne laser. 
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Fig. R— Contrast ratio and pulse degradation as functions of epsilon. 
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APPENDIX 



A general form for the expectation values of cos(x) such as the 
ones required for the calculation in the text is derived below. In par- 
ticular, it is shown by induction that 

<cos(x) \n) = n^J^il ( 19) 

where: x = ^"_j U , 

each t t is uniformly distributed on an interval [— e f , c ( ], the density 
function 

/,(/) - 1/(2*), /£[-e,-,e,.] 

= otherwise 

all of the U's are statistically independent. 

(cos (x) | n) is the expectation value of cos (x). 
We calculate the n = 1 case: 

(cos (x) | 1) = 7^- f ' cos (x) dx = ^iii). (20) 

2e,- J- ti Ci 

We now assume the form (19) and show that it can be extended to the 
n + 1 case: 

(cos (x) | n + 1) = / cos (x)F n+1 (x) dx 

J- oo 

where F„ +1 (x) is the distribution function of x for n + 1 terms. But 
F n+1 (x) = F n (x) i n+1 (x), 

= ^- / F n (x - r) dr. 
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Thus 

(cos (x) | n + 1) 

1 r°° r fn+ ' 

= / cos (x) I F n (x — t) dr dx 

"*n+\ J-oa J-tn+% 

= — — / / j cos (x) cos (t) -f- sin (x) sin (t)}F„(x) dx dr 

but the second term vanishes because F n (x), being generated by the 
successive convolution of even functions, is even. The integral reduces 
to 

(cos (x) I n + 1) = / cos (t) dr / cos (x)F n (x) dx 

= sin (t n+ i) fysinJO 

which completes the proof. 



n t 
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